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Abstract 

We have computed the form factors for B ^ it and D — > K{Tr) semileptonic decays 
on the lattice by using full non-perturbative 0{a) improvement, in the quenched 
approximation. Our results are expressed in terms of few parameters which describe 
the Q'^-dependence and normalization of the form factors. 
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1 Introduction 



In order to extract the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \Vub\ from the 
experimentally measured semileptonic branching ratio, B{B^ —>■ TT^i^u) it is essential 
to have a good theoretical control over the corresponding decay amplitude. This is difficult 
to solve from ffist principles because of non-perturbative hadronic effects. 

The relevant matrix element for a generic heavy (H) -flight (P) decay of pseudoscalar 
mesons is parameterized as 

{P{p)\V,\Hip^)) = (^p^ + p - qH^L^^ F^(g^) + ^^~^^ .Fo(g^) , (1) 

where = Q'yfj,q is the vector current involving the heavy [Q) and light (g) quark fields, 
-Fo/+(g^) are the form factors which carry the dynamical information on non-perturbative 
QCD, and q = {pn — p) is the momentum transferred from the parent to the daughter 
meson. 

In the case of massless leptons (which is an excellent approximation for i = e, n), the 
differential decay rate reads 



with A(g^) = (g^ + — mp)^ — Am^^m^p, being the usual triangle function. To obtain 
the total width, one should integrate (Q) over the entire (large) physical region, < < 
{niH — TTipY, which requires the precise knowledge of the normalization (say -F+(0)) and 
the g^-dependence of the form factor. 

There have been many studies devoted to the computation of these form factors. Var- 
ious quark models have been employed (see e.g. ref. P]), which in many aspects helped 
phenomenological understanding of the heavy— flight transitions. More quantitative pre- 
dictions, however, can be only obtained with approaches which are based on ffist QCD 
principles: QCD sum rules [] and lattice QCD simulations 

In this paper we present new results on the form factors -F+/o(g^), as obtained from our 
computation of the matrix element (0) on the lattice. In order to reduce lattice artifacts, we 
work on a large lattice (24^ x 64) with a small spacing (a ~ 0.07 fm, or = 2.7(1) GeV), 
using the non-perturbatively 0{a) improved Wilson action and operators P], i.e. all lead- 
ing lattice discretization effects proportional to the lattice spacing are absent. Due to the 
limitations on computational resources, however, one is obliged to make several approxi- 
mations, of which the most important ones are the use of the quenched approximation and 
the inability to resolve hadrons heavier than 3 GeV. The latter implies that the final answer 
for — > TT transition has to be reached through an extrapolation to rriB, which introduces 
large systematic uncertainties (which are of order 15%). An alternative is to discretize 
some effective theory, such as NRQCD, and to compute (|1|) from the opposite side, that 
is by extrapolating from infinitely heavy decaying mesons and eventually including higher 

^ The most recent results for i? — s- tt transition, obtained by using the Ught cone QCD sum rules 
(LCSR), can be found in ref. [|. 
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order I/ttih corrections [0. We also mention that alternative strategies to treat B —>■ Triu 
decay on the lattice have been developed in refs. [|, ||. 

The form factors are conveniently expressed in terms of three parameters {ch, olh-, Ph), 
which define the normalization and the dependence and which follow simple scaling laws 
as functions of the heavy quark mass. In terms of these parameters we have: 



Ch (1 - an) 



g2) (1 - anq"^) 



f„tf ) = SSJI^ , (3) 

where = q'^/mj^,, and rriH* is the lightest meson that couples to the vector current in 
eq. (|l|). This form is convenient to be used in Monte Carlo simulations for the experimental 
analyses of semileptonic decays. The values of the parameters are given below, in eq. (^). 
The parameterization which we refer to as BK, has been introduced in ref. [0, and 
it encodes most of the known constraints on the form factors. This includes: 

i) the kinematical constraint, -F+(0) = -Fo(O); 

ii) the heavy quark (meson) scaling laws predicted by the heavy quark effective theory 
(HQET), ~ y/ruH and Fq ~ Xj^froH-, which are applicable in the zero recoil 
region (g ^ 0, «.e. ^ g^^^) [|ll|]; 

iii) the heavy quark scaling law predicted by the large energy effective theory (LEET) 
and explicitly realized in the LCSR framework 0, -F+/o ~ 1/"^//^? applicable in the 
large recoil region, g^ ^ O, [11] ; 



iv) the knowledge of the position of the first pole in the crossed channel (g^ = mj^,), 
determining the g^-dependence of F^{q'^) in the small recoil region; 

v) the symmetry relation between the two form factors Fq = {2Ep/mH)F+, which is 
valid when the energy released to the light meson is large [|l3l |. 

In the absence of the last constraint v), one additional parameter would appear in eq. (|^) [|T0[] . 
However, we verified that the relation, Fq = {2Ep/mH)F+, holds when applied to our data 
(at large recoils), within the statistical accuracy. Its validity has also been verified in the 
LCSR approach 

In the infinite quark mass limit, the quantities (cHy/mn, (1 — C(H)fnH, {Ph — 1)"T'h) 
should scale as a constant. As usual, the l/mn and l/mjj corrections can be estimated 
from the fit with our lattice data. This allows us to extrapolate to the B meson mass. 
This strategy is employed in this paper and we call it Method I. Another way to handle 

^This scaling law has been pointed out long ago by Chernyak and Zhitnitsky in ref. [l^ . 
■^Very recently, it has been shown in ref. that this relation also holds to a high accuracy when 
radiative corrections are included. 
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the extrapolations has been recently proposed by the UKQCD collaboration [jl^. In this 
case the assumptions are different, namely the extrapolations are performed first in the 
light meson at fixed q^, and then in the heavy one at a fixed value oi v ■ p where v is 
the four-velocity of the heavy meson. The dependence of the extrapolated form factors 
is described by the parameterization @. We also employed this strategy, which we call 
Method II. Both methods lead to fully consistent results. They are also in a good agreement 
with those of UKQCD p, ||]. 

In the form (|), and by using the two different methods, we obtain: 

Method I : Method II : 

F(0) = 0.26(5)(4) F(0) = 0.28(6)(5) 

cs = 0.42(13)(4) cb = 0.51(8)(1) 



as = 0.40(15)(9) ae = 0.45(17)+-°^ 



(4) 

13 



where the first errors are statistical and the second are systematic. The result for -F(O) is 
displayed for the reader's convenience, although it is clear that F{0) = cb{1 — as)- We 
note that the results of the two methods are also in good agreement with the LCSR results 



of ref. |T6 |. 



Knowing the normalization and the g^-dependence of -F+(g^), we can integrate eq. 
to get the decay width. We quote: 

^ r(E° ^7r+£P) = (7.0±2.9) ps"^ (5) 



IK 



ub 



12 



obtained by combining the results of the two methods and adding the statistical and 
systematic errors in the quadrature. From this result, by using the measured branching 
ratio B{B^ — > 7r~^iu) = (1.8 ± 0.6) • 10"*^ |jl[] and the average meson lifetime r^o = 
1.548(32) ps [g, we get 

iKfel = (4.1±1.1)-10-^ (6) 

In addition, we also study Z)° —>■ ir^iu and Z)° K^lv decays and compare the results 
with the experimental data. 

The remainder of this paper is organized as follows: in sec. ^ we give the information 
on the lattice setup and describe the extraction of the raw form factor data; in sec. |] we 
explain the details used in extrapolation procedures to get the physical results which are 
presented in sec. ^ for B Tciu and in sec. ^ for D — > K{t[)Iv\ we briefly conclude in 
sec. |6|. 
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2 Lattice Setup and Computation of the Matrix Ele- 
ment 



In this section we give the essential information about the lattice setup used in our simu- 
lation. We also sketch the strategy to compute the matrix element (|I|) improved at 0{a). 



csw 



m 



1.614 



f^heavy — 



0.1250; 0.1220; 0.1190; 0.1150 
0.1344; 0.1349; 0.1352 
2.7(1) GeV 
0.13585(2) 



7(0) 



Cv 

bv 



22 



^,12 



0.79 
-0.09 
1.40 



Table 1: Parameters used in this work: the Wilson hopping parameters corresponding to the 
heavy (kq) and hght (Kg) quark masses; csw provides the non-perturbative 0{a) improvement 

of the Wilson action; cy, by, Zy^ ensure that the vector current is free of 0{a) effects (see 
eq. (0)). For each parameter we quote the reference where the quantity has been computed non- 



perturbatively. The critical hopping parameter {kc 
fixed as explained in ref. 



and the inverse lattice spacing (a ) are 



The results presented in this paper are obtained from a simulation on a lattice of size 
24^x64, using a sample of 200 independent gauge field configurations generated, at (3 = 6.2, 
in the quenched approximation. The values of the Wilson hopping parameters correspond- 
ing to the heavy and light quarks are listed in tab. |I]. In the same table, we give the values 
of the other parameters necessary for the 0{a) improvement which we implement in this 
work (note that all the parameters are determined non-perturbatively). The statistical 
errors are estimated by using the standard jackknife procedure, by decimating each time 
5 configurations from the whole sample. In our previous papers ||2^ , |2^ , we discussed the 
mass spectrum and decay constants using the 0{a) improvement as obtained on a subset of 
100 configurations ^. The updated values of our heavy-light pseudoscalar {Ha) and vector 
(H^) meson masses in physical units, obtained after extrapolating to the average up — down 
quark mass (see ref. for details) are: 



Mh, = { 0.642(2), 0.739(2), 0.829(2), 0.942(3) } 



The values of the heavy-light decay constants were updated in refs. 
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niHa [Ge V] = 



{ 1.74(7), 2.01(8), 2.25(9), 2.56(10) } 



Mjy. = { 0.684(3), 0.775(3), 0.860(3), 0.968(4) } 
mj^.[GeV] = { 1.86(7), 2.10(8), 2.33(9), 2.63(11) } . (7) 

As in our previous publications, we denote the meson masses in lattice units by capital 
letters, whereas the small case letters are used for the same masses in physical units, e.g. 
Mh = rriHCL. 

To compute the matrix element (0) on the lattice with improved Wilson fermions, the 
appropriate definition of the vector current reads ETI 



= Z'^y^ {l + ^{amq + amg) j Q-^^q + acy ■ id^Qa^^q 

with cr^j, = (i/2)[7^, 7jy], and where the renormalization constants Zy^ and by, as well as 
the subtraction constant cy, are known quantities (in tab. |I] we also quote the reference 
in which a given quantity has been computed non-perturbatively) . However, since we deal 
with the heavy-light operator, one can suspect that higher order terms, oc (am)" {n > 2), 
may be important. Their cancellation at tree level can be achieved by multiplying Zy by 
the KLM factor which we modify to include 0{a) improvement ||24|| , 



I.e. 



1 + am 

where arrig^Q} = {l/ng^Q} - 1/kc,.)/2 and m = {rrig + mQ)/2. 

To access the matrix element (|l]) on the lattice, one computes the following three- and 
two-point correlation functions 



C^\t.,ty-qjH) = j dxdy e'^^^-y-P^--'\P{0)V,{y)H\x)) , 

Cf}{t) = [ dx e'P-^J{x,t)j\0)) , (10) 



where H = qi'j5Q and P = qij^q are the source operators of the heavy-light and light- 
light pseudoscalar mesons, respectively, and J denotes any of the two {H or P). On the 
temporal axis in Cj^ \ we fixed H{x) to tr^ = 27, and the improved renormalized vector 
current V^{y) defined in eq. (P) is free (0 < < 63). pn is the three-momentum given 
to the decaying meson, and q denotes the momentum transferred to the daughter meson 
{q = ph — p). In this work we consider the following 13 kinematical configurations: 

'(0,0,0); (1,0,0); (1,1,0)} (11) 



Ph = 


(0,0,0) 




Ph = 


(1,0,0) 


Q 


Ph = 


(1,1,0) 
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where each component is given in units of elementary momentum on the lattice (27r/La). 
As usual, to get a better statistical quality of the signal, we averaged over the equivalent 
configurations, i. e. the ones that can be obtained by applying the cubic lattice symmetry, 
as well as the parity and charge conjugation transformations. 

From the asymptotic behavior of the three-point correlation function 



Cf\t.^y, qjn) ""^"^ ^ e~^^^y ■ {P{ph - mAm{Vn)) ■ ^ e-^-(*^-^) , (12) 

it is obvious that the removal of the exponential factors can be achieved by considering the 
ratio 

Cpp{ty,PH-q) C}j'H{U,-ty,pH) 

where y/ Zh = |(0|g'i75Q|iJ) | and y/Zp = \{0\qi'y5q\P)\, are obtained from the fit with 
asymptotic behavior of the two-point correlation functions. 

Cfjit) ^ §-e--'^ . (14) 

When the operators in the ratio (|T3|) are sufficiently separated, one observes the stable 
signal (plateau), which is the desired matrix element: 

RM'^^^"" {P{p)\V,\Hipn)) . (15) 

The plateaus are typically found for 12 < ty < 15. Fitting R^ty) to a constant in that 
interval, we extract the matrix element (p. The other possibility is to express R^ty) in 
terms of form factors for each time-slice, t = ty, and then to fit each form factor to a 
constant in the common plateau interval. As an important cross check of this procedure. 



we employ the following two methods 



- Ratio method consists in forming the ratio (|TBp of the correlation functions computed 
on the lattice, and using Zp and obtained from the fit with the corresponding 
two-point correlation functions in the same time-interval as the one at which the 
plateau has been attained. 

- Analytic method refers to the procedure in which we replace the two-point correlation 
functions in eq. (|T3|) by their asymptotic expressions ( [I^ ) and Zp and Zh are ex- 
tracted from the separate study of the two-point correlation functions. In this case 
we also take into account the symmetry of Cj j (t) with respect to the time inversion, 
t (64 — t). In this work we use the latticized free boson dispersion relation: 

sinh^(^)=smh^(f)+sm^(f) . (16) 

which appears to be appropriate in describing our data p3|, ^ . 
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Figure 1: Plateaus of the two form factors for the specific combination of the heavy and light 
quark masses (corresponding to hq = 0.1220 and Kg = 0.1349). The form factors are extracted 
from the fit in t G [12,15]. 



The differences in the results obtained by using the ratio and analytic methods are always 
smaller than the statistical errors. In fig. 0, we illustrate the form factors Fo/+(t) as 
functions of time t = ty, obtained from the ratio (|13D by using both the ratio and the 
analytic methods. From now on, our central results will always be those obtained by using 
the ratio method. The difference between that result and the one we get by applying 
the analytic method, will be added to the systematic uncertainties. We have also checked 
whether the results change after replacing \p\ 2sin(|p|/2). This difference, which is 
an O(a^) effect, is around 1% and we neglect it. This is in agreement with findings of 
ref. where this source of discretization errors was studied in great details for the 
gluon propagator and the triple gluon vertex. These errors, although important for large 
momentum injections, were shown to be completely negligible for the small momenta (such 

This indication makes us confident that this 



as the ones considered in this paper (11 



source of systematic errors is under control. 



3 Getting to the physical results: Extrapolations 



Once we computed the form factors for all 13 momentum combinations (11), we need to 



address the question of the mass extrapolations. Firstly, since our light meson masses are 
in the range 0.54 GeV < rrip < 0.84 GeV, a small extrapolation to the physical kaon 
and a longer one to the pion mass are needed. Secondly, the extrapolation in the heavy 
quark mass is needed, if we are to reach the physical results relevant for B ^ n decay. 
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-0.2 0.0 0.2 0.4 



Figure 2: The form factors Fo/+(g^) for Hq = 0.1349 and hq = 0.1220, plotted as functions of 
/My, where My is the mass of the hghtest vector meson exchanged in the i-channel. The hnes 
show the fit to the parameterization (|3|). 



These extrapolations are closely related to the fit of the form factors because both 
extrapolations, in light and heavy mesons, obviously enlarge the physical region accessible 
from the semileptonic decay, G [0, (m/f — mp)^]. Before discussing the extrapolations, we 
list in tab. ^the BK-parameters (|^) resulting from the fit to our data for each combination 
of the heavy and light quarks. One such a fit is shown in fig. |^. 

3.1 Method I: Extrapolation of the parameters 

The parameters (pi = {ch, an, Ph), for each pair of the heavy and light hopping parameters 
(see tab. ^ give us information on the normalization and the g^-dependences of both form 
factors. For fixed heavy quark, one can apply the lattice plane method to extrapolate 



to the case when the final meson is either a kaon or a pion. This essentially means a linear 
fit in the light quark mass (Mp oc rriq) for each parameter: 

iPi = A + Bi Ml , (17) 

and then an extrapolation to M^, determined in the (Mp,Mi/) plane by fixing the ratio 
Mv/Mp to the physical point mp/m.,^. Similarly, the value of M^- is determined from 
the physical rnK'/friK ratio. The result of that extrapolation is illustrated in fig. 0, and 
the extrapolated parameters are listed in tab. |[ The form factor at = is just the 
combination c/f(l — an)- We note that whether we extrapolate the parameters ch and 
an to the final pion (kaon) separately and then combine them to F{0), or we combine the 
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Figure 3: Extrapolation of each of the parameters {ch, an, Ph) to the final kaon/pion. The 
plot is made for the heavy quark corresponding to kq = 0.1220. 
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Ch 


an 






izoU — 


io44 


L.zoyLZ ) 






U.o i[Z) 


1250- 


1349 


1.09(13) 


0.29(9) 


1.37(12) 


0.77(3) 


1250 - 


1352 


0.99(13) 


0.29(11) 


1.33(13) 


0.70(4) 


izzU — 


1 Q /I /I 

io44 




n oo /'^^ 
U.Zol^y j 


1 /io/'i Q^ 
i.4z(^io j 




1220 - 


1349 


0.99(12) 


0.26(10) 


1.41(13) 


0.73(3) 


1220 - 


1352 


0.88(11) 


0.26(12) 


1.35(14) 


0.65(4) 


T T on 


T Q /I /I 

io44 




U.Zol^y j 


i.4o(^io j 


U. r 0(^0 j 


1190- 


1349 


0.93(11) 


0.26(11) 


1.42(14) 


0.69(3) 


1190 - 


1352 


0.83(10) 


0.27(11) 


1.33(12) 


0.61(4) 


1150 - 


1344 


1.02(11) 


0.28(10) 


1.43(14) 


0.73(3) 


1150- 


1349 


0.88(10) 


0.26(12) 


1.42(14) 


0.65(4) 


1150 - 


1352 


0.77(9) 


0.25(12) 


1.35(13) 


0.58(4) 



Table 2: The results of the fit of our form factor data to the parameterization given in eq. @, 
for all combinations of nq-Kg. 

parameters to obtain F{0) and then extrapolate to the pion (kaon), the result remains the 
same. The latter number is also given in tab. ^. 

After the extrapolation in the light mass, one needs to fit the parameters in the inverse 
heavy meson mass and to interpolate to the D or extrapolate to the B meson mass. The 
quantities which scale in the infinitely heavy quark limit (0), are fitted as 

0e jci^^W^, (l-aH)Mj,, (/5h-1)Mh} , (18) 

where we neglect the logarithmic corrections. In fig. §we show the 1/Mh behavior of each 
parameter for the case in which the final light meson is the pion. Again, as a cross check, 
one can also fit F{Q)M^^'^ to the form (0). The results for -F(O) extrapolated to Mb in 
that way is fully consistent with -F(O) obtained by combining separately the extrapolated 
values of cb and to -F(O) = cb(1 — a^). As can be seen from fig. ^, however, a systematic 
difference in the extrapolated values of the parameters is obtained depending on whether 
we perform a linear (62 = 0) or a quadratic (62 free) fit in the inverse heavy meson mass. 
Indeed, besides the quenched approximation, this difference represents the main systematic 
uncertainty present in our calculation. To make a comparison with other approaches easier, 
we give the values of the l/rriH corrections (in physical units) that are obtained from the 
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0.0 I — ^ — ' — ' — ' — ' — ' — ' — ' — ' — ' — ^ 




Figure 4: Extrapolation in the heavy meson mass. The curves represent the results of the 
quadratic (full) and linear (dashed) extrapolations. 



12 



Final meson — kaon Final meson — pion 

, s , ^ s 











F(0) 








F(0) 


0.1250 


0.96(14) 


0.29(11) 


1.33(13) 


0.69(4) 


0.82(16) 


0.29(14) 


1.30(17) 


0.59(6) 


0.1220 


0.85(12) 


0.25(12) 


1.36(14) 


0.64(4) 


0.71(12) 


0.24(15) 


1.32(17) 


0.55(6) 


0.1190 


0.80(11) 


0.26(12) 


1.34(13) 


0.60(4) 


0.66(11) 


0.25(14) 


1.29(15) 


0.51(6) 


0.1150 


0.74(10) 


0.24(13) 


1.35(13) 


0.57(4) 


0.60(10) 


0.23(15) 


1.31(15) 


0.48(6) 



Table 3: Values of the BK-parameters after extrapolating to the final vr and K meson. The 
extrapolation of the form factor at = provides a cross check of the extrapolation. 



linear and quadratic fit with our data: 

^, , (3.1 ±0.5) GeV=^/2 
F(0) = ^ X 



m 



3/2 
H 



(0.98 ±0.09) GeV 



f (0) = '^-^^l-l'"-^''^ X 



m 



3/2 
H 



1 - 



mu 

(2.0 ±0.3) GeV 

mu 



± 



;i.2±0.1 GeV)^ 



m 



H 



.(19) 



We observe, in passing, that the l/m^ corrections are large and comparable to the ones 



which appear with the heavy-light decay constants |24 



At this point we should decide which result, from either the linear or quadratic fit, 
to quote as the central one. As a criterion for that, we test whether the result of the 
extrapolation satisfies the soft pion theorem, i.e. the Callan-Treiman relation applied to 
the case of 5 ^ tt^z/ 0: Q 



B 



(20) 



where fs and /tt are the corresponding decay constants. The values of the ratio fs/ fn 
will be given in the next section. It turns out that after the quadratic extrapolation in 
1/Mh (|11), we find a better consistency of Fo(m^) with fs/ f-K- Thus, the result of the 
quadratic l/M/^-extrapolation will be taken as our central value, and the difference between 
the quadratically and linearly extrapolated values will be accounted for in the estimate of 
the systematic error. All the numerical results of the extrapolations will be given in next 
section, where we discuss our physical results. 



3.2 Method II: UKQCD 

The second method that we use in this paper is the one proposed by the UKQCD col- 



laboration in ref. [17|. As in the previous subsection, the starting point are the results 



^An extensive high statistic study of this relation on the lattice, in the NRQCD framework, has been 
made in ref. 0. 
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given in tab. ^ For a given value of the heavy quark mass, one extrapolates the form 
factors to the final pion (kaon) by using (pT]). In order to separate the intrinsic light quark 
mass dependence of the form factors from the one that arises from the variation in g^, the 



extrapolation in the light quark mass should be made at fixed q [17|. In other words, we 



consider the extrapolation to the pion for a constant v ■ p, which is defined as 

' ^''^ 

where v is the four- velocity of the heavy meson and p is the momentum of the light meson. 
We note that the lattice results for the form factors are in the region 0.2 < v ■ p < 0.5. 
In this window we could fix 7 equidistant points v ■ p = 0.20, 0.25, . . . 0.50, where the 
distance between the points is taken to be the same as in ref. ||T^. We used the fit to the 
parameterization Q to interpolate to the chosen values of v ■ p. For each fixed v ■ p, 
we then extrapolate to the final pion as described in the previous subsection. A typical 
situation is displayed in fig. ^, where the heavy quark corresponds to kq = 0.1220, i.e. 
Mh^ = 0.739(2) in lattice units. The result of the extrapolation is marked by the filled 
circles in fig. |^. 

We now have to extrapolate in the heavy quark mass for which one uses the heavy 
quark scaling laws. For each fixed value of f ■ p, we extrapolate 



4(M„) = + ^ + ^ . (22) 



where 



The term multiplying both form factors gives the logarithmic dependence on the heavy 
meson mass predicted by the HQET in which the vector current has the anomalous di- 
mension whose leading order coefficient is 70 = —4 |3^ 0. As usual, the extrapolation in 
l/Mfj ( p2|) is performed either linearly (02 = 0), or quadratically (02 free). Obviously, the 
HQET scaling law [|T^] is valid in the region close to the zero- recoil {v ■ p —>■ 0), and the 
extrapolation from the points which are rather far from zero-recoil is an assumption. The 
difference between the results obtained after the linear and quadratic extrapolations will 
be included in the systematic uncertainty. In fig. we show this extrapolation for the 
case of f ■ p = 0.35. As in the previous subsection we quote the result of the quadratic 
extrapolation as our central value, and the difference with the linearly extrapolated value is 
included in the estimate of the systematic error. The resulting form factors for the B tt 
semileptonic decay are given in tab. ^. The last step is the fit of the data from tab. ^ to 
the form (|^), which is what we discuss next. 

^Note that the hght meson mass in eq. is the physical pion mass. 

^Other interpolation forms are also applicable, such as the pole formula (with the pole mass left as a 
free parameter), which is what has been done in ref. [0. We checked that the final results on the form 
factors do not depend on the form used in this interpolation. 

^It should be noted, however, that the effect of the logarithmic corrections in eq. (|2^) on the extrapo- 
lation to ms is practically negligible. 
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Figure 5: Lattice data for the form factors Fq/_^_{v-p) with the heavy quark fixed to hq = 0.1220, 
and for three different values of light quark masses. Dotted (full) curves represent the fit to the 
parameterization (^). Filled symbols denote the form factors extrapolated to the final pion at a 
given value of v ■ p which has been fixed according to eq. (^). 
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0.5 1.0 1.5 0.5 1.0 1.5 

Figure 6: Extrapolation of tiie form factors ( p2[ ) at v ■ p = 0.35 in the inverse heavy meson mass 
to the Bd meson. Dashed and full curves correspond to the linear and quadratic extrapolation, 
respectively. 



g2 [GeV'] 






13.6 
15.0 
16.4 
17.9 
19.3 
20.7 
22.1 


0.46(7)+:°^ 
0.49(7)^:°^ 
0.54(6)1:°^ 
0.59(6)^:?^ 
0.64(6)^:?^ 
0.71(6)1;?^ 
0.80(6)1:?^ 


0.70(9)^:J° 

0.79(10)t:J2 
0.90(10)+;^° 

1.05(11)1:^° 

1.25(13)+;°^ 
1.53(17)^;°^ 
1.96(23)1;°^ 



Table 4: The values of the semileptonic B ^ it form factors at several values of . First errors 
are the statistical and the second are systematic whose estimate is explained in the text. 
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Quantity 


Method 1 Method U LCbR 16 


Cb 

CtB 

Pb 


0.42(13)(4) 0.51(8)(1) 0.41(12) 
0.40(15)(9) 0.45(17)1;?^ 0.32tfr 
1.22(14)1:^2 1.20(13)1;^^ - 




0.26(5)(4) 0.28(6)(5) 0.28(5) 

1.3(6)1:° 1-5(5)1:° - 
1.21(9)1:°° 1-19(11)1:?? 1-281:}^ 


\Vui,\^^-V{B^ ^-K+lv) [ps-i] 
IK^I/IO-^ 


6.3 ±2.4 ±1.6 8.5 ±3.8 ±2.2 7.3 ± 2.5 
4.3 ±0.9 ±0.6 4.0 ±1.0 ±0.5 4.0 ±0.7 ±0.7 



Table 5: Main results of this paper, as obtained by using the two methods explained in the text. 
For easier comparison, we also list the results obtained by using the LCSR |16| ]. 



4 Physical results 

4.1 B ^ iriui 

The main physical results as obtained by employing the Method I and Method II, are listed 
in tab. ^. In the same table we also give the results obtained by using the LCSR which 



are fitted to the parameterization (0) in ref. [16]. First, we see that the final results of the 



two methods agree very well with each other. We also notice a very good agreement with 
the results reported in ref. [|T^ where the Method 71 has been used. This is shown in fig. |^. 
Moreover, we see that there is a pleasant agreement with the LCSR results. We verified that 
the Callan-Treiman relation ( pOD is well satisfied when the quadratic extrapolation in the 
heavy meson mass is performed. We recall that on the same set of configurations, we have 
/b/Ztt = 1.29(10) and 1.46(16), as obtained from the linear, and quadratic extrapolation 
of the ratio fn/fn in '^/Mh, respectively. The latter value has to be compared with the 
value of Fo(m^) given in tab. ^ 

We were also able to examine the effect of the SU(3) breaking at = 0, by replacing 
the final pion by a kaon. This ratio (see tab. ^ turns out to be large compared to one and, 
although with large errors, quite compatible with the findings of the LCSR. 

We now discuss the sources of systematic errors which are combined in the quadrature. 

- Extraction of the matrix element: The central results are obtained by using the ratio 
method (see sec. The whole analysis has also been performed with the matrix el- 
ement extracted by using the analytic method. The difference in the results is always 
smaller than the statistical error, and we include it in the systematic uncertainty. 

- Discretization errors: The main source of discretization errors is expected to be of 
(9((amj|/)2), where mn is the heavy meson mass. Our central results are obtained 
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Figure 7: The form factors relevant for B Tiiv decay obtained by using Method II, are directly 
compared to the result of the UKQCD group (only statistical errors are shown). We also draw 
the curves describing the dependence with the parameters given in tab. |. 
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by correcting the 0{a) improved renormalization constant by the KLM factor as 
given in eq. (^. We repeated the analysis without the KLM factor and found only 
a small discrepancy with the central results of the parameter cb (of about 5 %) 
which is included in the systematics. The other parameters remain unchanged. We 
emphasize that this comparison can only give us idea about the size of discretization 
errors. To check this point we computed Zy for the case of the heavy-heavy vector 
current Q'joQ for each of our heavy quarks. By considering the ratio 



dxdy{H{y) Q7oQ(x,t) H^{0)) 

n{t) = , (24) 

2 dx {H{x,t)W{0)) 



which is shown in fig. H, one directly accedes l/Zy- We follow ref. [17], and denote 



such an obtained Zy as Zy^^^ which we now compare to the values that one gets by 
using the expression Zy = Zy\l + byarriQ), where Zy^ and by are those given in 
tab. I: 



t^Q = 


0.1250; 




= 1.155(31) [Zy 


= 1.144(6)], 


I^Q = 


0.1220; 




= 1.254(35) [Zy 


= 1.252(6)], 


I^Q = 


0.1190; 


^y 


= 1.363(40) [Zy 


= 1.367(6)] , 


I^Q = 


0.1150; 


^y 


= 1.516(47) [Zy -- 


= 1.528(6)]. 



(25) 

Clearly, results of the two procedures are fully consistent with each other. In other 
words, there is no room for a term proportional to (amg)^. 

Besides, we also used the value of a~^(//^) = 2.8(1) GeV and repeated the whole 
analysis which leads to the results only slightly different (less than 5%) from our 
central values. This difference is, however, added in the systematic uncertainty. 

- Heavy meson extrapolation: As already mentioned, we quote as central values the re- 
sults obtained by using the quadratic extrapolation in 1/Mh. The difference between 
this and the results of the linear extrapolation is added to the systematics. Beside 
quenching, this is the main source of systematic uncertainty present in our analysis 
(it is of order 15%). 

The small differences of BK-parameters given in tab. ^give a more pronounced effect in 
the spectrum of the differential decay rate (0). This is shown in fig. After integrating 
eq. (0) we get the decay width (divided by jV^bj^), which is also given in tab. |^. As 
mentioned in the introduction, after comparing this result to the measured branching ratio 

^From the fit Zy'^^^ — Zy\l + byaniQ + by^^\amQY), where the coefScients Zy^ and by are given 
in tab. |l|, we obtain by^^^ — 0.05(27), thus consistent with zero. 
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Figure 8: I/Z'^ff^ as obtained from the ratio 7^(^) (see eq. (p^) for the time slices between the 
two source operators (placed at and 27 respectively). Ratio is plotted for each heavy quark 
mass 2mQ = 1/kq — 1/K.cr- Fit is performed for t E [5,25]. 



o amQ=0.320 
□ amQ=0.418 
• amQ=0.521 
■ amQ=0.667 
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B{B^ — > -R^lv) |jT|, we are able to make an estimate of |V^b| f^. We also checked that 
if the pole/dipole form is used to describe the g^-dependence of the form factors that we 
obtained by using Method II, the result is fully consistent with the one we quote in tab. H, 
namely we obtain \Vub\~^ ■ F(i?° n+iiy) = 7.8 ± 3.2 ± 2.1 ps^^ 

To make a unique estimate, we combined in quadrature the statistical and systematic 
errors for F(5° —>■ ir'^iu), and made a weighted average of the results of two 

methods to get, |Kt;,|^^ F(i?° — > n^iu) = (7.0 ± 2.9) ps. Now, as before, if we compare 
this result to the measured branching ratio, we obtain \Vub\ = (4.1 ± 1.1) • 10^'^, which is 
consistent with the world averaged value from inclusive decays, presented at this year's 
ICHEP conference pi: |Kfe| = (4.1 ± 0.6 ± 0.2) ■ lO'^. 



5 D ^ TTiui and D Klu^ 

In this section we briefly discuss the results relevant for the semileptonic decays of D meson. 
As in the previous section, we apply both Method I and Method II. In this case, however, 
the most important source of systematic uncertainty is negligible, namely in order to reach 
the physical value of the heavy meson mass, only a smooth interpolation of the data from 
tab. 1^, by using eq. ([T^), is required. Notice that for these decays, a bulk of final lattice 

The last number, however, should be taken cautiously because the experimental result has been 
obtained by relying on various different theoretic estimates. 
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Figure 9: Differential decay rate for B nil' decay. We sliow the curves as obtained by both 



methods. We also plot the result of ref. |16], obtained using the LCSR. 



points is out of the physical region (this is illustrated in fig. |TOD. The values of the form 
factors whose g^'s are within the physical region for the semileptonic decays are listed in 
tab. §. 

Our final results are given in tab. where we again compare to the numbers obtained 
from the LCSR study of ref. [TB|. 

As in the previous section, we plot the dependence of the partial decay width for 
both decays {D ^ K and D ^ tc). The agreement with the results of the LCSR is less 
good than it was for B TTiu. Our central numbers for the form factors at = are 
smaller. This region (of small g^'s) dominates when integrating over the phase space. On 
the other hand the discrepancy in the g^-dependence of the form factors is less important 
since the kinematically accessible region from these decays is small. Our results for the 
parameter a, larger than those obtained by LCSR, suggest that the contribution of the 
higher excited states to the form factor -F+(g^) is not absent. LCSR, instead, predict that 
the g^-dependence of F^^''^^\q^) is very close to a pure pole behavior (with the pole mass 
being m^* ). We also mention that cd gives access to the value of the strong coupling 
of D and tt to the first pole D* (known as gD*D-K)i namely Qd'O-k = "^cditld* / Jd* . With 
fn* = 258(14)(6) MeV |2§, this amounts to gD-Dn = 12 ± 2 and go-Dn = 14 ± 5, by 
using the Method I and Method II, respectively. The most recent LCSR prediction is 
go'Dn = 10.5 ± 3 10 0. A detailed discussion on this coupling with a complete list of 



The quoted LCSR results for D — > K£iy Q are obtained by using ms{fJ, = rric) = 150 MeV. 
^^The coupling go'D-n is one of the main parameters in the approach based on the phenomenological 



lagrangians |35 
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Figure 10: Form factors for D ^ K and D ^ it semileptonic decays. The full line correspond 
to the results of Method I, and the dashed one to the results of Method II. The physical region 
< q'^ < (rriD — ^k/tt)'^ is indicated in both cases. 



g2 [GeV^] 


-0.04 


0.47 


0.97 


1.48 


F^--{q') 


0.56(6)^:°2 
0.56(6)1:°^ 


0.62(6)1:°2 
0.67(6)1:°J 


0.70(6)1;°S 
0.81(7)^::°^ 


0.80(6)1:°J 
1.03(9)1:°J 




0.19 


0.69 


1.20 


1.70 


F^-^iq') 


0.68(4)(0) 
0.70(5)(0) 


0.76(4)(0) 
0.84(5)(0) 


0.85(4)(0) 
1.02(6)(0) 


0.96(4)(0) 
1.29(7)(0) 



Table 6: The semileptonic D ^ it and D ^ K form factors at several g^'s obtained by using 
Method II. Only the values at belonging to the physical region for the semileptonic decays are 
listed, with the exception of = —0.04 GeV^ which is very close to = 0. The systematic 
errors are estimated in the way explained in the text. 
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Figure 11: The same as in fig. (^) but for D meson semileptonic decays. 
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Quantity 


Method 1 Method U LCbR [16 


cdk 

OLDK 


0.90(12) 1.17(23) — 
0.27(11)1:°? 0.43(12)1:°° -0.07tfr 

1.34(14)1:°^ 1.22(13)1:°^ - 


CDtt 


0.76(13)1:°^ 0.88(18)1:°° 
0.27(14)1:°° 0.36(16)1:°° o.oii:J^ 
1.31(17)1:^^ 1.23(17)1:^^ - 




0.66(4)1:°J 0.66(4)1:°J 0.78(11) 
0.57(6)1:°J 0.57(6)1:°^ 0.65(11) 

1.6(3)1:° 1-9(5)1:° 

1.16(5)1:[!5^ 1.22(6) ~ 1.2 



Table 7: The same as in tab. |. 



references can be found in refs. | p6|| . 

Finally, in tab. ^ we also give our results for the total decay rates, which are in a fair 
agreement with the experimental values |]TU|. As in the previous subsection, we can also 



1-2 



IT lU 



compare our result for | Vcs \ 
K-iu) = (8.5 ± 0.5) ■ 10-2ps-i. We obtain 



to the experimental measurement T^'^^'^'\D'^ 



IK, 



1.07 ±0.09 



(26) 



which is to be compared to \Vcs\ = 1.04±0.16 [19|, extracted from the D ^ K semileptonic 



decay where the vector meson dominance for the form factor F+(g^) has been assumed. 

m 



the 



Both central values are larger than 1, although consistent with \Vcs\ = 0.9745(8 
value obtained by using the unitarity of the CKM matrix. The more precise experimental 
result for B{D° K 
clarifying the situation. 



?z/) expected from CLEO and E791 |Q will hopefully help in 



6 Conclusion 

In this paper we have computed the form factors relevant for the semileptonic decays of the 
heavy-to-light pseudoscalar mesons on the lattice. We used the fully 0{a) improved Wilson 
action and vector current. Since the actual computation is performed on a lattice with an 
inverse lattice spacing smaller than 3 GeV, the physical results for the most interesting 
decay in this class, namely B iriu, had to be reached through an extrapolation in the 
inverse heavy meson mass, for which HQET and LEET provide us helpful scaling laws. To 
describe the interplay of the mass dependence and the g^-dependence of the form factors. 
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i [L) — >■ K lu) [ps J 


i [U — > 71 IV) [ps J 


Experiment ||T9[ 


(8.5 ±0.5) ■ 10-2 


(8.6 ±1.8) -10-3 


This work {Method I) 
This work {Method IT) 


(6.9 ± 1.1) ■ 10-2 
(7.3 ± 1.1) ■ 10-2 


(5.6 ± 1.4) ■ 10-3 
(5.9 ±1.5) -10-3 


APE 95 [0 
UKQCD 95 H 
LANL 95 


(9.1 ±2.0) ■ 10-2 
(7.0 ± 1.7) ■ 10-2 
(9.4 ± 1.5) ■ 10-2 


(8±2) ■ 10-3 
(5.2 ±1.8)- 10-3 
(7.4 ± 1.1) ■ 10-3 


LCSR |16 


(9.2 ±3.5) ■ 10-2 


(6.3 ±2.4) ■ 10-3 



Table 8: Comparison of the fuh decay widths obtained in this work, to the experimental value, 
to the several recent lattice results and to the prediction of ref. |jl6[ . In the computation we used 
\Vcs\ = 0.9745(8) and \Vcd\ = 0.220(3) |l9|. 



we used two different methods which are described in the text. We summarize our findings 
as follows: 

- There is a good agreement of the results obtained by using two different methods to 
reach the physically interesting results of our study; 



We note a very good agreement with the results of ref. [0 in which the same (im- 
proved) lattice action has been used; 

Our results agree with the ones of ref. |jl6| where the LCSR were employed. This is 



especially the case for B — > vr^z/ decay. We also note that, although with somewhat 
large errors, we see SU(3) breaking effects for the form factors at g2 = Q, consistent 
with the findings of the LCSR studies; 

- We also verified that our form factor data are consistent with the soft pion theorem 
(Callan-Treiman relation (^0|) ). 



The main results of our study are given in the introduction and in tables A remaining 
uncertainty for the B n decay arises from our inability to work with the B meson 
directly on the lattice. This uncertainty can be further reduced by repeating the calculation 
on a larger and finer lattice. Besides, our lattice simulation has been performed in the 
quenched approximation. It is obviously desirable to attempt the unquenched study of the 
semileptonic decays. 
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